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ABSTRACT 



o 
in 



The most promising astrophysical sources of kHz gravitational waves (GWs) are 
the inspiral and merger of binary neutron star(NS) /black hole systems. Maximizing 

■ the scientific return of a GW detection will require identifying a coincident electro- 
J> ' magnetic (EM) counterpart. One of the most likely sources of isotropic EM emission 

ON , from compact object mergers is a supernova-like transient powered by the radioactive 

decay of heavy elements synthesized in ejecta from the merger. We present the first cal- 
culations of the optical transients from compact object mergers that self-consistently 
determine the radioactive heating by means of a nuclear reaction network; using this 
heating rate, we model the light curve with a one dimensional Monte Carlo radiation 
<0 transfer calculation. For an ejecta mass ~ 10 _2 Mq [10~ 3 -Mq] the resulting light curve 

peaks on a timescale ~ 1 day at a V-band luminosity vL v ~ 3 x 10 41 [10 41 ] ergs s _1 
(My = — 15[— 14]); this corresponds to an effective "f" parameter ~ 3 x 10 -6 in the 
Li-Paczynski toy model. We argue that these results are relatively insensitive to un- 
certainties in the relevant nuclear physics and to the precise early-time dynamics and 
ejecta composition. Since NS merger transients peak at a luminosity that is a factor 

■ ~ 10 3 higher than a typical nova, we propose naming these events "kilo-novae." Due 

to the rapid evolution and low luminosity of NS merger transients, EM counterpart 
searches triggered by GW detections will require close collaboration between the GW 
and astronomical communities. NS merger transients may also be detectable following 
a short-duration Gamma-Ray Burst or "blindly" with present or upcoming optical 
transient surveys. Because the emission produced by NS merger ejecta is powered by 
the formation of rare r-process elements, current optical transient surveys can directly 
constrain the unknown origin of the heaviest elements in the Universe. 

Key words: binaries: close; supernovae: general; stars: neutron; gamma rays: bursts; 
gravitation; nuclear reactions, nucleosynthesis, abundances 



1 INTRODUCTION 

The direct detection of gravitational waves (GWs) would 
be a major breakthrough for both fundamental physics 
* E-mail: bmetzger@astro.princeton.edu and astrophysics. With upgrades of the ground-based in- 
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terfe rometers LIGQ (Abramovici et al.l ll992T ) and Virgo 
(e.g. ICaron et al] 1 19991 )" to "advanced" sensitivity expected 
within the next decade, GW detection is rapidly becoming 
a realistic— even anticipated— possibility. 

The most promising astrophysical sources of GWs for 
ground-based detectors are thought to be the GW-driven 
in-spiral and coalescence of binary compact objects (neu- 
tron stars [NSs] and black holes [BHs]). Advanc es in gen- 
eral r elativistic simulations of the merger process l|Pretoriusl 
120051 ) are honing our understandi ng of the strength and 
form of the expected signal (see e.g. lFabedl2009l ; |Puej[2009l . 
for recent reviews). However, estimates of the merger rates 
based on known NS-NS binaries and population synthe- 
sis remain uncertain by at least an ord e r of magnitude 
jKim et al] 120051; iBelczvnski et al] 120061; iKalogera et al ' 



2007; LIGQ Scientific Co llaboratio n fc Virgo Collaboration 
20101 ) For instance. iKim et al] (120051 ) estimate that the NS- 



NS merger rate detectable with advanced LIGO will be 
27^21 P er y ear 5 implying that if the true rate lies on the 
low end of present estimates then only a few sources may be 
detected per year. This possibility makes it especially crucial 
that we extract the most science from each event. 

Optimizing the science from a detected GW signal re- 
quires ident ifying a coincident electromagnetic (EM) coun- 
terpart (e.g.|Schutzlll986l.l2002l;ISvlvest rc 2003; S tubbs!l2008l ; 
iBloom et alj|2009al ; |Phinnevll2009l ; IStamatikos et alj|2009h . 
By independently identifying the source's position and time, 
several of the degeneracies associa t ed with the GW sig- 
nal are lifted l|Hughes fc Holzll2003l ; lArun et al] 120091 ) and 
the sign al-to-noise required for a c onfident detection is de- 
creased (|Kochanek fe Pirad fl993l ; iDalal et~alll2006l ). Cou- 
pled with its GW-measured luminosity distance, identify- 
ing the merger's redshift (e.g. by localizing its host galaxy) 
would also allo w for a precision measurement of the Hubble 
constant (e.g. iKrolak fc Schutzlll987l ; iHolz fc Hughesll2005l ; 
iDeffavet fc Menoull20o'7i r Likewise, the potential wealth of 
complementary information encoded in the EM signal may 
be essenti al to fully unra veling the astrophysical context of 
the event (|Phinnevll2009l ). 

The most commonly discussed EM signal associated 
with NS-NS/NS-BH mergers is a short-duration Gamma- 
Ray Burst (GRB), powered by the accretion of material 
that remains in a centrifug ally-supported torus around the 
BH f ollowing the merger JPaczvnskil Il98d; iNaravan et al] 
1 19921 ). The Swift satellite ([Gelrrels et all |2004| ) has rev- 
olutionized our understanding of short GRBs by detect- 
ing and localizing a significant number of their afterglows 
for the first time. This has enabled the discovery that 
short GRBs likely originate from a more evolved stel- 
lar population than those of long-duration GRBs (e.g . 
IBloom et alj|2006l ; lBerger et al] |2005: Vill asenor et al . 2005; 



Bergeij|2009f for a recent review), consistent with an ori- 



gin a ssociated with compact object mergers |Nakar et al] 
1 2006! ) ■ Despite these suggestive hints, however, it is not 
yet established that all short GRBs a re uniquely associ- 
ated with NS-NS/NS-BH mergers (e.g . iHurlev et al1l2005l ; 
iMetzger. Quataert fc Thompson] [2008) nor that all merg- 
ers lead to an energetic GRB. Furthermore , only a smal l 
fraction of GRB jets are pointed towards us |Rhoadslll999l ) 
and for off-axis events, the prompt and afterglow emission 
are much dimmer due to relativistic de-beaming. Although 
some emission may be observed by off-axis viewers, such "or- 



phan" afterglows (e.g. iTotani fc Panaitescu] |2002) are typi- 
cally expected to peak at radio wavelengths on a timescale 
of months-years (ISoderberg et al.l 120061 ; iRossi et al] 120081 ; 
IZhang fc MacFadven 20091 ). Only a limited fraction of short 
GRBs are detected in radio, even when viewed on-axis (e.g. 
ISoderberg et al]|2006l ). 

In parallel to the advances in GW detectors, the ad- 
vent of large-scale optical surveys with increasing sensitiv- 
ity, rapid cadence, and sky-area coverage is leading to a 
revolution in the study of transient obje cts. These i nclud e 
the Palomar Transient Factory (PTF; iRau et al.l I2009TI . 
Pan- STARRs faaiser et al]|2002l ). SkyMappe r dKeller et al 



2007), and the VLT Survey Telescope (VST; iMancini et al 
2000J) , which are paving the way for future end eavors such as 



the L arge Synoptic Survey Telescope (LSST; Istrauss et al.l 
120101 ) and the pr oposed Synoptic All Sky Infrared Imaging 
(SASIR) survey (|Bloom et al.ll2009bl ). Given these present 
and anticipated future capabilities, the most promising 
EM counterpart of compact object mergers is arguably an 
isotropic, optical/near infrared (NIR) wavelength signal. In 
addition to providing time-stamped maps of the night sky 
for use in follow-up observations, these "blind" surveys could 
also detect EM counterparts even independent of a GW or 
GRB trigger (see ^3|) . 

One proposed source of relatively isotropic opti- 
cal/NIR emission following a NS-NS/NS-BH merger is a 
supernova(SN)-like transi ent powered by the radioactive 
decay of merger ejecta dLi fc Paczvriskil 1 19981; hereafte r 
LP98; cf. lKulkarnill2005l ; IMetzger. Piro fc Quataert] |2008| ). 
Although Type la superno va light curves are power ed largely 
by the decay of 56 Ni (e.g. iKasen fc Wooslevll2009l ). most of 
the ejecta from NS-NS/NS-BH mergers is highly neutron- 
rich (electron fraction Y e ~ 0.1 — 0.4) and produces little 
Ni. Instead, much heavier radioactive elements are formed 
via rapid neutron capture (r-process) nucleosynthesis fol- 
lowing the decompression of the ejecta from nuclear den- 
sities (e.g. lLattimer fc Schramml 1 1974 Il97rj : lEichler et all 
ll989l : lFreiburghaus et al.lll999l ). Although the r— process i - 
self lasts only a matter of seconds, these newly-synthesized 
elements undergo nuclear fission, alpha and beta decays on 
much longer timescales as they descend to /3-stability. The 
resulting energy release can power detectable thermal emis- 
sion once the ejecta expands sufficiently that photons can es- 
cape. Due to the lower quantity of ejecta and its faster speed, 
however, the resulting transient is dimmer and evolves faster 
than a normal SN. Transients from NS mergers are thus a 
challenge to detect and identify. 

Although the basic LP98 model provides a qualitative 
picture of the thermal transients from NS-NS/NS-BH merg- 
ers, it makes a number of simplifying assumption and leaves 
several free parameters unconstrained, including the fraction 
of nuclear energy released and the precise distribution of 
decaying nuclei. LP98 further assume that the photosphere 
radiates as a black body, which is a poor assumption at mod- 
erate optical depths and in light of the substantial UV line 
blanketing that may accompany the rich energy spectra of 
the very heavy nuclei that dominate the composition. These 
details may be important for predicting the unique, "smok- 
ing gun" features of merger-related transients. Because the 
transient sky is e xpected to be rich in its diversity (e.g. 
iBecker et afll2004l ). more detailed predictions may be essen- 
tial to identifying candidate sources in real-time for deeper 
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follow-up observations, especially considering the likelihood 
that only limited information (e.g. photometric colors) may 
be available. Understanding the detailed spectroscopic prop- 
erties of merger transients is clearly an important endeavor. 

In this work we present the first self-consistent calcu- 
lations of the optical/NIR counterparts to NS-NS/NS-BH 
mergers. In particular, our work goes beyond previous work 
in two important ways: (1) we use a nuclear physics reaction 
network to calculate the radioactive heating of the ejecta and 
(2) we employ the Monte Carlo radiative transfer code SE- 
DONA to more accurately model the light curve and colors 
of the resulting EM transient. We begin in |J2]with prelim- 
inary considerations, including a discussion of the sources 
of neutron-rich ejecta from NS-NS/NS-BH mergers f ffTTj) 
and a brief review of the relevant scalings for radioactively- 
powered transients ( i)2.2[) . In §3] we describe the nucleosyn- 
thesis that occurs as the material decompresses from nu- 
clear densities and our calculations of the resulting radioac- 
tive heating, including a detailed discussion of the efficiency 
of fission//?— decay thermalization f M3.2|) . In §31 we present 
calculations of the light curves and color evolution of NS- 
NS/NS-BH merger transients, highlighting the unique fea- 
tures of these events and the primary uncertainties in the 
theoretical predictions. We find that the peak luminosities 
of NS merger transients are typically ~ few xlO 41 ergs s _1 , 
or a factor ~ 10 3 larger than the Eddington luminosity for 
a solar mass object. We therefore dub these events "kilo- 
novae," since standard novae are approximately Eddington- 
limited events. In §5] we discuss the implications of our re- 
sults for the present constraints on, and the future detection 
of, kilonovae from NS-NS/NS-BH mergers, including the di- 
rect constraints that optical transient surveys place on the 
astrophysical origin of r-process elements f i|5.3.1[) . We sum- 
marize our results and conclude in §6] 



2 PRELIMINARY CONSIDERATIONS 

2.1 Sources of Neutron- Rich Ejecta 

There are several potential sources of neutron-rich ejecta 
from NS-NS/NS-BH mergers. First, neutron-rich material 
can be ejected due to tidal forces during the me r ger itself 
( Lattimer fe Schramm! Il974l ; iRosswog et alj 1 19991 ; Rosswod 
l2005h . The quantity of this dynamically ejected material de- 
pends sensitively on the NS-NS/NS- BH binary par ameters 
and the NS equation of state (e.g. iRosswoel l2005h . Since 
this material primarily originates from the NS's neutron-rich 
oute r core, it has a typical elec tr on fraction V e ~ 0.03 — 0.1 
(e.g. lHaensel fc Zdunikl [l990al lb1; iRosswod l2005h . The elec- 
tron fraction probably remains low since the ejecta remains 
cold (and hence thermal weak interactions remain slow) due 
to adiabatic losses as the material rapidly expands from nu- 
clear densities (e.g. iRuffert et al.lli"997l ; iDuez et alj|2009h . A 
typical outflow speed is v ~ 0.1 c. 

Neutron-rich material also originates from outflows 
from the accretion disk on longer, viscous timescales. 
Neutrino-heated winds are driven from the disk for 
a variety of accretion rates and disk radii during 
its early evolution (e.g. iMetzger. Piro fc Quataerj 120081 ; 
ISurman et al. 1 120081 ; iDessart et al.l 120091 ). Although these 
outflows are generally ne utron-rich, they can be proton- 



2008; iBarzilav fc Levinsonl l2008'). An even larger quantity 
of mass loss occurs at later times once neutrino cool- 
ing shuts off, due to powerful outflows driven by viscous 
heating and the nuclear recombination of free nuclei into 
a— particles dMetzger. Piro Sz Quataert l200Sl: Metzger et al.l 
l2009al ; iLee et al l l2009h . IMetzger et all l|2009al ) show that 
~ 20 — 50% of the initial disk mass is ejected with a range 
of electron fractions Y e ~ 0.1 — 0.4. The wind's asymptotic 
speed in this case is also v ~ 0.1 — 0.2 c, set by the ~ 8 MeV 
per nucleon released as heavy elements are formed. 

In summary, considering both the tidally- and wind- 
driven ejecta from NS-NS/NS-BH mergers, from M e j ~ 
up to ~ O.IMq in neutron-rich ejecta is expected with 
v ~ 0.1 c and Y e <C 0.2, i.e. sufficiently neutron-rich to un- 
dergo a robust (low entrop y) third-peak r-proces s during its 
subsequent expansion (e.g. Hoff man et al.lll997h . A similar 
amount of material may be ejected with Y e ~ 0.2 — 0.4. 
Although this material is not sufficiently neutron-rich to 
reach the third r-process peak, it also produces heavy ele- 
ments that contribute a comparable radioactive heating rate 
(Fig.©. 



2.2 Analytic Estimates 

The majority of the energy releas ed by the r-process occurs 
on a timescale of ~ seconds (e.g. Frciburg haus et all 1 19991 ; 
iGorielv et al.ll2005l ; IMetzger et alJl201Ch . However, most of 
this initial heating (and any residual heat from the merger 
itself) is lost to adiabatic expansion because the outflow is 
highly optically thick at these early times. A significant EM 
luminosity is only possible once the density decreases suffi- 
ciently th at photons ca n escape the ejecta on the expansion 
timescale (|Arnettlll982T l. The photon diffusion time through 
the outflow at radius R is approximately 



t d = 



BnM e 



cR 



(1) 



wher e k is the opacity an d B ~ 0.07 for a spherical outflow 
(e.g. iPadmanabhanl [2000 V Setting this equal to the expan- 
sion time i cxp = R/v gives the characteristic radius for the 
EM emission to peak 



R 



( BvkM c 



peak 



1/2 



1.2 x 10 cm 



(o.lc) 



1/2 



W- 2 M Q 



1/2 



(2) 



where we have taken k — 0.1 cm g _1 as an estimate of the 
line opacity of the r-proce ss ejecta, assuming it is similar to 
that of Fe-peak elements (|Pinto fc EastmanTl2000r i . We dis- 



cuss the validity of this assumption further in §471] Assum- 
ing free expa nsion R = vt, Rpe^ is reached on a timescale 
|Arnettlll982h 



ipeak ^ 0.5 



days (ay 



-1/2 



io- 2 m 



1/2 



(3) 



rich in some cases (e.g. IMetzger. Thompson fc Quataerj 



The above expression is strictly valid only if t pea k exceeds 
the intrinsic radioactive decay lifetime of the ejecta. This 
condition is generally satisfied for r— process ejecta due to 
their rather short /3-decay half-lives. This short timescale 
tpoak ~ 1 day compared to that of a normal SN (t poa k ~ 
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Figure 1. Radioactive heating rate per unit mass E in NS 
merger ejecta due to the decay of r-process material, calculated 
for the Y e = 0.1 ejecta trajectory from Rosswog et al. (1999) 
and Freiburghaus et al. (1999). The total heating rate is shown 
with a solid line and is divided into contributions from j3— decays 
(dotted line) and fission (dashed line). For comparison we also 
show the heating rate per unit mass produced by the decay chain 
56 Ni 56 Co 56 Fe (dot-dashed line). Note that on the ~ day 



timescales of interest for merger transients (t 



caki 



ED 



fission and /3— decays make similar contributions to the total re- 
process heating, and that the r-process and 56 Ni heating rates 
are similar. 



weeks) is one of the defining characteristics of kilonovae from 
NS mergers. 

Provided that the radioactive power can be approxi- 
mated as a decreasing power-law in time Q oc t~ a with 
q < 2, the brightness of the event depends most sensitively 
on the amount of radioactive heating that occurs around 
the timescale t pea k: Q pea k = L Qdt w Q(t P cak)i P eak = 



L Qdt 

Jt pcak 

/M e jC 2 , where / < 1 is a dimensionless number (LP98). 
Parametrized thus, the peak bolometric luminosity is ap- 
proximately 



^pcak 



/peak 



td(Rpc 



""""-""■'(^(iffi) 



1/2 



1O" 2 M 



1/2 



(4) 



and the effective temperature is given by 



/ r \ ^ 4 

/ -^peak 



1.4 x 



-1/8 



1O- 2 M 



-1/8 



(•») 



Note that L pea k oc /, yet the value of / is left as a free 
parameter in the LP98 model, with values up to / ~ 10~ 3 
considered plausible a priori. In £|3.2l we present explicit cal- 
culations of Q and show that the effective value of / is 
~ 3 x 10~ 6 . Thus, for M c j ~ 10 _2 Mq we expect a transient 
with peak luminosity ~ 10 42 ergs s _1 (bolometric magni- 
tude M bo i ~ —16) and a photospheric temperature ~ 10 4 K, 
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Figure 2. Final abundance distribution from the fiducial model 
with Y e = 0.1 (Fig. [TJ, shown as the mass fraction versus 
mass number A. Measured solar system r-process abundances are 
shown for comparison with black dots. They are arbitrarily nor- 
malized to the computed abundances for A = 195. 



corresponding to a spectral peak at optical/near-UV wave- 
lengths. 



3 RADIOACTIVE HEATING 
3.1 Network Calculations 

In this section we present calculations of the radioactive 
heating of the ejecta. We use a dynamical r-pr ocess network 
l|Martinez-Pinedcll2008llPetermann et all2008h that includes 
neutron captures, photodissociations, /3— decays, a— decays 
and fission reactions. The latter includes contributions from 
neutron induced fission, /3 delayed fission, and spontaneous 
fission. The neutron captur e rates for nuclei with Z ^ 83 are 
obtained from the work of iRauscher fc Thielemann! (|200fj) 
and are based on two differe nt nuclear mass m odels: the 
Finite Range Droplet Model l|Moller et all Il995l ) and the 
Quenched version of the Extended Thomas Fermi wit h 
Strutinsky Integral model (ETFSI-Q) (|Pearson et al.lll996t ). 
For nuclei with Z > 83 the neutron captu re rates and 
neutro n-induced fission rates are obtained from lPanov et all 
(2009). Beta-decay rates includi ng emission o f up to 3 neu- 
trons after beta decay are from iMoller et all (|2003D . Beta- 
delayed fission and sp ontaneous fission rates are deter- 
mined as explained by iMartmez-Pinedo et all ((2007) . Ex- 
perimental rates for alpha and beta decay have been ob- 
tained from the NUDAT database^ Fission yields for all 
fission proces ses are determined using the statistical code 
ABLA l|Gaimard fc Schmidtll99ll ; lBenlliure et alll998T ). All 
heating is self-consistently added to the entropy of the fluid 
following the procedure of Freiburghaus et al. (1999). The 
change of temperature during the initial expansion is de- 
termined using the Timmes equation of state (Timmes & 
Arnett 1999), which is valid below the density p ~ 3 x 10 11 
g cm -3 at which our calculation begins. 

As in the r-process calculations performed by 
Freiburghaus et al. (1999), we use a Lagrangian density p(t) 



1 http : //www.nndc .bnl .gov/nudat2/ 
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Figure 3. Total radioactive heating rate per unit mass E, cal- 
culated for several values of the electron fraction Y e of the ejecta 
and for different nuclear mass m odels (see text). Ca lculations em- 
ploying the FRDM mass model l|Moller et aljl995l ) are shown for 
Y e = 0.05 (triple dot-dashed line), Y e = 0.1 (solid line), Y e = 0.2 
(dotted line), and Y e = 0.3 (dot-dashed lin e). A calculation em- 
ploying the EFTSI-Q jPearson et al.lll996l) mass model is shown 
for Y e = 0.2 (dashed line). Note that on timescales of hours— days, 
the radioactive heating rates in all models agrees to within a fac- 
tor ~ 4. 

taken from the NS-NS merger simulations of Rosswog et 
al. (1999). In addition to p(t), the initial temperature T, 
electron fraction Y e , and seed nuclei properties (A,Z) are 
specified for a given calculation. We assume an initial tem- 
perature T = 6 x 10 9 K, although the subsequent r-process 
heating is not particularly sensitive to this choice because 
any initial thermal energy is rapidly lost to PdV work dur- 
ing the initial expansion before the r-process begins (Meyer 
1989; Freiburghaus et al. 1999). For our fiducial model we 
also assume Y e — 0.1, Z ~ 36, A ~ 118 (e.g. Freiburghaus 
et al. 1999). 

Our results for the total radioactive power E with time 
are shown in Figure [1] On timescales of interest the radioac- 
tive power can be divided into two contributions: fission and 
/3— decays, which are denoted by dashed and dotted lines, re- 
spectively. The large heating rate at very early times is due 
to the r-process, which ends when neutrons are exhausted 
at t ~ 1 s ~ 10~ 5 days. The heating on longer timescales 
results from the synthesized isotopes decaying back to sta- 
bility. On the timescales of interest for powering EM emis- 
sion (tp 0a k ~ hours-days; eq. [3]), most of the fission results 
from the spontaneous fission of nuclei with A ~ 230 — 280. 
This releases energy in the form of the kinetic energy of 
the daughter nuclei and fast neutrons, with a modest con- 
tribution from gamma-rays. The other source of radioactive 
heating is fi— decays of r-process product nuclei and fission 
daughters (see Table [T] for examples corresponding to our 
fiducial model) . In Figure [1] we also show for comparison 
the radioactive power resulting from an identical mass of 
56 Ni and its daughter 56 Co. Note that (coincidentally) the 
radioactive power of the r-process ejecta and 56 Ni/ 56 Co arc 
comparable on timescales ~ 1 day. 

In Figure [5] we show the final abundance distribution 
from our fiducial model, which shows the expected strong 
second and third r-process peaks at A ~ 130 and A ~ 195, 



respectively. For comparison, we show the measured solar 
system r-process abundances with points. The computed 
abundances are ra t her d ifferent to the one obtained by 
iFreiburghaus et alj (| 1999h due to an improved treatment of 
fission yields and freeze-out effects. 

Although we assume Y e = 0.1 in our fiducial model, 
the ejecta from NS mergers will possess a range of electron 
fractions (see i)2.ip . To explore the sensitivity of our results 
to the ejecta composition we have run identical calculations 
of the radioactive heating, but varying the electron fraction 
in the range Y e = 0.05 — 0.35. Although in reality portions of 
the ejecta with different compositions will undergo different 
expansion histories, in order to make a direct comparison 
we use the same density trajectory p(t) as was described 
earlier for the Y e — 0.1 case. Figure [3] shows the heating 
rate for ejecta with Y e = 0.05, 0.2 and 0.3 in comparison to 
the fiducial model with Y e = 0.1. Although the heating rate 
for different values of Y e differs substantially at early times 
(< 10~ 4 days), E agrees between the models to better than 
a factor ~ 2 at the later times that are the most important 
for transient EM emission. 

Our results for E could in principle also be sensitive to 
the assumed properties of the nuclei in the r-process path 
(e.g. masses and neutron-capture cross sections), which are 
uncertain and must be obtained via theoretical modeling. 
In our fiducial model (Fig. [TJ we employ the FRDM model 
(Moller et al. 1995) for nuclear masses. In order to explore 
the sensitivity of our results to the assumed nuclear physics, 
we also performed an oth erwise identical calc ulation using 
the ETFSI-Q mass model dPearson et al.|[l996l ). as shown in 
Figure [3] for Y e = 0.2. Although the two models again differ 
in their early-time predictions for E, on timescales > 1 hour 
they converge to a heating rate within a factor < 4. 

Finally, although the Lagrangian density trajectory p(t) 
that we employ in our fiducial model formally corresponds 
to dynamically-ejected rather than wind-driven ejecta, both 
are likely present in NS-NS/NS-BH mergers (see fl2~Tj) . 
Thus, we have also performed an otherwise identical cal- 
culation, but instead using a trajectory p(t) appropriate for 
(higher entrop y ) disk win ds, similar to those studied in e.g . 
Arcones et ail (120071) (cf. iMetzger. Piro fc Quataertl 120081 ; 
Surman et al. I l2008l ). Although we do not show our results 
explicitly in this case, we find that the heating rate E de- 
creases in a similar manner to the dynamically-ejected ma- 
terial on timescales ~ 1 day. However, the overall normaliza- 
tion of E is smaller by a factor ~ 10 because in high entropy 
winds the mass fraction of heavy nuclei is reduced at the ex- 
pense of a higher alph a particle fraction, w hich do not con- 
tribute to the heating (|Hoffman et al.|[l997h . Although some 
of the wind-driven material in NS mergers may have high en- 
tropy i|Surman et alJl2006tlMetzger et al.ll2008D . most of the 
total mass ej ected likely has low-m odest entropy (S ~ 3 — 10 
kb barvon~ 1 ; lMetzger et al. 2009a). When correcting our re- 
sults for the higher mass fraction of heavy nuclei in a lower 
entropy wind, we find that the heating rate E on timescales 
~ 1 day in the wind ejecta agrees within a factor ~ 2 to 
that of the dynamically-ejected material. 

To summarize, the heating rate for our fiducial model 
in Figure [T] (which we employ throughout the remainder of 
the paper) appears to be relatively insensitive to the precise 
trajectory and composition of the ejecta, and to the uncer- 
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tainties in the nuclear properties of the unstable nuclei near 
the r-process path. 

In order to understand why we find such a robust heat- 
ing rate on timescales ~ 1 day, it is first instructive to com- 
pare r-process ejecta with that produced in Type la SNe. 
In Type la SNe, the ejected material is processed through 
nuclear statistical equilibrium with Y e ~ 0.5. This favours 
the production of N — Z nuclei and, in particular, 56 Ni. The 
56 Ni nucleus (N = Z = 28) is produced in high abundance 
both because 28 is a magic nucleon number and because 
even-even (N = Z) nuclei have an additional binding energy, 
commonly known as the "Wigner energy." At the high tem- 
peratures at which 56 Ni is produced, atoms are fully ionized 
and, consequentially, 56 Ni cannot decay by atomic electron 
capture. In this case th e half-life has been computed to be 
ty 2 ~ 4 x 10 4 years by iFisker et all |T999). Once the tem- 
perature decreases sufficiently that the inner _R"-shell orbit 
electrons recombine, the decay proc eeds at the laborato ry 
measured rate t 1 / 2 ^ 6.075(10) days (|da Cruz et al.lfl99"3 h 

The situation is different for neutron-rich r-process 
ejecta. First, r-process nuclei decay by /3~ and hence the 
half-life is unaffected by the ionization state of the mat- 
ter. Secondly, the r-process results in a rather broad dis- 
tribution of nuclei with mass number spanning the range 
A ~ 110 — 210. Because the nuclei produced in NS merg- 
ers likely follow a distribution similar to their solar system 
abundances (see Fig. [2]), maxima will occur at the second 
(A ~ 130) and third (A ~ 195) r-process peaks. The overall 
r-process abundances peak in our calculations (as in the so- 
lar system) near the second peak, which is why second-peak 
nuclei dominate the j3— decay heating rate (see Table [1]). 

We argue below, however, that the energy generation 
rate E is approximately independent of the precise distribu- 
tion of heavy nuclei, provided that the heating is not domi- 
nated by a few decay chains and that statistical arguments 
can be applied. This conclusion is supported by our results 
in Figure [3] which show that E is relatively insensitive to 
the composition of the ejecta, despite the fact that differ- 
ent electron fractions can result in rather different abun- 
dance distributions. Perhaps most striking, the heating rate 
is similar whether the second r-process peak is produced 
via the fission of nuclei near the magic neutron numbers 
TV = 184 (A ~ 280), as occurs for highly neutron-rich ejecta 
(Ye < 0.2), or whether it is produced directly with little or 
no fission cycling, as occurs for Y e ;> 0.3. 

Assuming a broad distribution of exponentially decay- 
ing nuclei with mass number A the evolution of the energy 
generation rate E can be understood by the following argu- 
ments. For an isotopic chain of odd-A nuclei the Q-values 
are essentially proportional to the neutron excess r\ = N — Z 
and the beta decay rate A oc rf due to the 3-body nature 
of the final state. The situation is slightly more complicated 
for even-A chains due to the presence of pairing that in- 
creases the binding energy of even-even nuclei and modifies 
the global proportionality of the Q- value and neutron excess. 
However, the selection rules of beta decay favor a maximum 
change in angular momentum between initial and final states 
of one unit; as a result, the decay of odd-odd nuclei, that typ- 
ically have angular momentum J > 1, proceeds via excited 
states in the daughter even-even nucleus. Consequently, the 
global dependence A oc rf is recovered. Assuming that the 
number of nuclei per neutron excess interval is constant the 



Table 1. Properties of the dominant /3— decay nuclei at t ~ 1 day 



Isotope t 1/2 Q(«0 4 b) el C) 4 d) £" Vg<0) 
(h) (MeV) (MeV) 



135j 


6.57 


2.65 


0.18 


0.18 


0.64 


1.17 


129 Sb 


4.4 


2.38 


0.22 


0.22 


0.55 


0.86 


128 Sb 


9.0 


4.39 


0.14 


0.14 


0.73 


0.66 


129 Te 


1.16 


1.47 


0.48 


0.48 


0.04 


0.22 


132j 


2.30 


3.58 


0.19 


0.19 


0.62 


0.77 


135 Xe 


9.14 


1.15 


0.38 


0.40 


0.22 


0.26 


127 Sn 


2.1 


3.2 


0.24 


0.23 


0.53 


0.92 


134j 


0.88 


1.2 


0.20 


0.19 


0.61 


0.86 


56 Ni (/) 


146 


2.14 


0.10 


0.10 


0.80 


0.53 



(a) Total energy released in the decay; (b),(c),(d) Fraction of the 
decay energy released in electrons, neutrinos, and 7— rays; (e) 
Average photon energy produced in the decay; (f) Note: 56 Ni is 
not produced by the r-process and is only shown for comparison 
(although a small abundance of 56 Ni may be produced in accre- 
tion disk outflows from NS-NS/NS-BH mergers; Metzger, Piro, & 
Quataert 2008). 

number of nuclei per decay rate interval, A, and lifetime in- 
terval, r = 1/A, are given by: 

dN = N \~ i/B d\, dN = N T- 6/5 dT, (6) 

where No is a normalization constant. Further assuming that 
the energy generation at a time t is dominated by nuclei with 
t — t that release energy Q oc r _1//5 , the energy generation 
rate then becomes: 

E oc t~ 7/5 = t~ 1A . (7) 

The same result can be obtained assuming that we have 
a distribution of nuclei that follows equation , each releas- 
ing an energy Q oc A 1 / 5 with a rate A. In this case the energy 
generation rate is: 
/> 00 

Eoc XX 1/5 e - xt N \" i/5 dX = 7V r(7/5)r 7/5 , (8) 
Jo 

where F is the Gamma function. 

The above discussion neglects the fact that with increas- 
ing neutron excess the beta-decay populates an increasing 
number of states in the daughter nucleus. Consequently, we 
expect an exponent slightly larger than 5 for the depen- 
dence of decay rates with neutron excess. This will result 
in a power-law decay with an exponent smaller than the 
value of 1.4 deduced above. Overall, this analytic deriva- 
tion is in reasonable agreement with the numerical results 
in Figure[T] which correspond to £ oc t~ a with a ~ 1.1— 1.3 
on timescales of hours— days. Incidentally, we note that this 
functional form is remarkably similar to the heating rate 
E oc t~ 1,2 fou nd for the decay of nuclear waste from terres- 
trial reactors (|Cottingham fc Greenwoodll200ll : pg. 126). 

3.2 Thermalizing Processes 

Of the total power released by nuclear reactions E (Figure 
[1} , only a fraction ethcrm will thermalize with the plasma and 
hence be useful for powering EM emission. In this section 
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we estimate e t hcrm- Since the light curve peaks on approxi- 
mately the timescale i pea k (eq. [3]), we shall normalize our 
considerations to this time. 



3.2.1 /3-decay Heating 

First, consider the energy released by /3— decays, which dom- 
inate E at late times (Fig. [TJ. The total energy released in 
the decay Q is divided between the outgoing neutrino and 
electron, and the gamma-rays produced as the daughter nu- 
cleus cascades to the ground state from excited nuclear lev- 
els. In Table [JJ we list the properties of a sample of nuclei 
which contribute appreciably to E at t ~ t poa k ~ 1 day, 
as determined from our network calculations in £|3. 1 1 The 
properties listed include the decay half- life ti/2, the rela- 
tive fraction of Q carried away by the electron, neutrino, 
and gamma-rays (e e , e„, and e 7 , respectively) and the mean 
gamma-ray energy _E 7 vg . Most of this information was ob- 
tained or calculated using data from the Lawrence Berkeley 
National Laboratory's Isotopes ProjectQ 

Because the energy imparted to the outgoing electron 
E e = e e Q is generally greater than or similar to the electron 
rest mass (0.511 MeV), the electron is mildly relativistic 
and, as a result, carries a similar fraction of the outgoing 
energy as the neutrino (i.e. e e ~ e v ~). Although the neutrino 
readily escapes the ejecta and does not contribute to the 
heating, the electron is charged and interacts electromagnet- 
ically with the ambient electrons and nuclei. The dominant 
thermalizing process is electron-electron coulomb scattering. 

For electron-electron scattering in the fast test particle 
limit (E e S> kT, where T ~ 10 4 K is the temperature of the 
background plasma), the energy exchange (or "thermaliza- 
tion") timescale is given by 

C c e rm a 4.6 x 10 13 s (— f^) hiA-V" 1 , (9) 

where In A is the Cou lomb logarithm (e.g. NRL Plasma For- 
mularly; iHubal 12007). Assuming a spherical, homogeneous 
outflow, the electron density n e at time t is approximately 
given by 



e (47T/3)i?V 

where jj, e » Am n /Z is the mean mass per electron, m n is the 
mass of a nucleon, and we have assumed an average charge 
Z ~ 60 and mass A ~ 130 for the r-process nuclei. Although 
the r-process nuclei are only partially ionized on timescales 
t ~ tpcak, n e includes both free and bound electrons because, 
for purposes of high energy scattering, they have identical 
cross sections (the impact parameter for E e ~ MeV is much 
smaller than the atomic scale). 

Thus, the ratio of the thermalization time due to 
electron-electron collisions t^ e to the timescale at which 
the emission peaks is given by 

^thermApcak ~ 

2 http://ie.lbl.gov/education/ 




Equation (|lll) shows that for a typical value E e ~ 0.5 MeV, 
*therm > Vak for t <; 10 t pcak , implying that the /3-decay 
electrons will efficiently thermalize on the timescales of in- 
terest. 

Table Q] shows that typically ~ 1/2 of the /3— decay 
energy is also released in the form of ~ MeV gamma- 
rays. Although a portion of the gamma-rays will thermalize 
via Compton scattering and photoelec tric absorption (e.g. 
I Colgate et all Il980l ; ISwartz et all 119951 ), a significant frac- 
tion will also escape, especially as the ejecta expands and 
the optical d epth decrease s . In t he case of 56 Ni and 56 Co, 
for example, ISwartz et alj (|l995l ) find an effective absorp- 
tive opacity which is about ~ 15 per cent of the fully- 
ionized Thomson opacity (i.e. k 7 rss 0.03 cm 2 g _1 ) using 
Monte Carlo radiative transfer calculations. Since t pca k is 
attained when the optical depth is r pca k ~ c/v ~ 10, the 
thermalization "optical depth" is < 1 for times greater than 
« (K/K 7 )(u/c)t pcak ~ t pcak (see eqs. [T] and [3]). As Table 
[JJ illustrates, the mean 7— ray energy E* vg from 56 Ni de- 
cay is within a factor ~ 2 of those produced by the other 
j3— decays, so we expect similar 7— ray thermalization prop- 
erties in Type la SNe and in NS merger ejecta. We conclude 
that photons will partially thermalize for t < i pca k, but they 
will contribute little heating at later times (t 3> i pea k). 

From Table [TJ we infer average values of e e ~ e„ « 0.25 
and e 7 ~ 0.5 (Table [TJ. Combining our results, we conclude 
that the effective /3— decay thermalization fraction will vary 
from ethcrm ~ e e + e 7 m 0.75 to ethe™ ~ £ e ~ 0.25 as the 
ejecta expands from R <C -R poa k to R> i? poa k- 



3.2.2 Fission Heating 

In the case of fission most of the radioactive energy is re- 
leased as kinetic energy of the fission product nuclei, with 
a typical daughter energy of Ea ~ 100 MeV. In this case 
the dominant thermalizing process is coulomb scattering off 
ambient nuclei of similar mass A and charge Z. I n the case o f 
ion-ion collisions the thermlization timescale is !|Huball2007f ) 

t A ~ A ~ 

therm ^ 

5 x 10- s taA-v ( ea f 2 f av /2 f £\r 2) 

A VlOOMeV/ V 130/ Veo/ * 

where ua ~ p/Am n is the number density of ambient nu- 
clei. Thus, using equations (J3]) and (|10|) the ratio of the ther- 
malization timescale to the timescale at which the emission 
peaks is given by 

W m /Wk-6xl0 (^Mev) [ m ) X 

(irmur(i) 3 (i3) 

Since t ti ^ lm -C t pca k we conclude that the fission daughters 
will also thermalize on timescales ~ t pea k, implying that 
ethcrm ~ 1 for fission. Therefore, even though fission con- 
tributes less to E than /3— decays at t ~ t pca k ~ 1 day, its 
higher thermalized fraction suggests that it may dominate 
the heating. 
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3.2.3 Neutron Heating 

Neutrons, emitted either spontaneously following fission or 
induced by j3— decays, also carry a modest portion of the 
released energy. Although the contribution of neutrons to E 
is generally much less than that of j3— decays or the kinetic 
energy of fission daughters, we consider their thermalization 
as well for completeness. 

At high densities (e.g. in terrestrial reactors) fission- 
product neutrons can be captured by heavy nuclei to induce 
further reactions. For the much lower densities of present in- 
terest, however, the neutron capture timescale (with a typ- 
ical cross section ~ millibarns) is much longer than the free 
neutron beta-decay timescale ~ f 5 minutes. As a result, fast 
neutrons created by fission rapidly decay into into fast pro- 
tons (with a typical energy E p ~ 1 MeV) before capturing. 
In order to thermalize, the protons must exchange energy 
with the much heavier ambient charged particles (the pro- 
ton density is much too low for p-p scattering to be efficient). 
We find that the proton's thermalization time is larger by a 
factor ~ Z 2 A 1/2 (E p /E A f /2 ~ 10 2 than that of the fission 
daughter nuclei (eq. Q2]). This suggests that the proton will 
have ttherm ~ ipcak and hence may not thermalize. 

Our estimates above neglect, however, the possible ef- 
fects of magnetic fields, which can trap charged particles 
and enhance their thermalization if th e field is directed per - 
pendicular to the outflow velocity (e.g. IColgate et ahl ff98C0. 
For instance, if the NS involved in the merger has a (mod- 
est) surface field strength of B ~ fO 9 G, this translates into 
a field strength of B ~ 3 x 10 -8 G at R ~ i? pca k by flux 
freezing. The larmor radius for a 1 MeV proton at i? pca k is 
tl ~ fO 12 cm, which is <C -R pca k ~ 10 14 cm. This suggests 
that the proton's residence time (and hence thermalization) 
may be significantly enhanced due to the magnetic field. As 
a result, we conclude that the energy released in neutrons 
will also likely thermalize. 



3.2.4 Net Heating Efficiency and the Effective Value of f 

Considering both f3— decays (with e t hcrm ~ 0.25 — 0.75) and 
fission (ethcrm ~ 1), we conclude that the net heating frac- 
tion of merger ejecta is between e t hcrm ~ 0.25 and « 1, de- 
pending on time and the relative contributions of /3— decays 
and fission to E. 

From Figure [T] we find that E decreases approximately 
as a power law E cx t~ a , with a value a ~ 1.1 — 1.4 on 
timescales of hours— days, relatively close to the heating 
functional form adopted by LP98: Q cx fc 2 /t (see eq. [3] 
and surrounding discussion) . In £|3 . 1 1 we presented a sim- 
ple derivation of this result that explains why a <J 1.4 
(see eqs. [7]- ED- For an average thermalization efficiency 
of etherm ~ 0.75, our results imply that the effective value of 
/ is ~ 3 x 10 -6 at t ~ 1 day. This is somewhat lower than 
the range of values considered by LP98 and much lower than 
has been estima ted elsewhere in the literature. For instance, 
iRosswod J 2005) estimates a NS merger transient peak lu- 
minosity Lpcak ~ 10 44 ergs s _1 , corresponding to an ef- 
fective value of / - 10~ 3 for M cj ~ 10~ 2 A/ Q (eq. g]). 
He derives this under the assumption that an appreciable 
fraction of the total energy released in forming heavy r- 
process nuclei (B/A ~ 8 MeV nucleon -1 ) is released over 
a timescale t pca k ~ 1 day. This is incorrect because most 



of the binding energy is released by the formation of seed 
nuclei in the initial expansion (on timescales ~ millisec- 
onds) and by the subsequent r-process (on timescales < 1 
second; see Fig. [TJ. Since this heating all occurs at radii 
R £ (v x 1 second) ~ (f/0.3c)10 10 cm, this early-time heat- 
ing suffers a factor > 10 4 loss in thermal energy due to PdV 
work before the outflow expands to the radius i? pca k ~ 10 14 
cm (eq. [2]) at which photons can finally escape. Instead, 
the luminosity at t ~ t pca k is primarily powered by residual 
energy released as r-process products fission and decay back 
to stability. This occurs on much longer timescales and in- 
volves a significantly smaller energy release (typically closer 
to ~ 10 -3 MeV nucleon -1 on timescales ~ 1 day). 



4 LIGHT CURVES 

4.1 Radiative Transfer Calculation 

We calculate the light curves of compact object merger 
ejecta u sing the time-depen dent radiative transfer code SE- 
DONA (|Kasen et all 120061 ). Our set-up is similar to that 
described in Darbha et al. (2010, in prep) for the case of 
56 Ni decay-powered transients produced by the accretion- 
induced collapse of white dwarfs. However, in the present 
case we have modified the code to include the radioactive 
heating Q = e t hcrm-E for NS mergers as calculated in ij3] 
Although SEDONA can track 7— ray thermalization, we do 
not use this option given the large number of decaying nu- 
clei and their complex 7-ray spectra; rather, we incorporate 
the escape fraction into an average thermalization efficiency 
ethcrm that we hold constant in time. We calculate mod- 
els assuming both ethcrm = 0.5 and ethcrm = 1 in order to 
bracket the uncertainty in the precise fraction of 7— rays that 
thermalize (see H3.2.4p . 

Although the ejecta from NS-NS/NS-BH mergers is 
likely to be highly asymmetric (e.g. the "banana-like" ge- 
ometry of tidal tails; Rosswog 2005), we assume a one- 
dimensional (spherically symmetric) geometry for simplicity. 
We shall address the effects of the full multi-dimensional 
kinematics of the outflow in future work. Since SEDONA 
uses a velocity-time grid, a required input to the calculation 
is the velocity profile p(v) of the (assumed homologous) ex- 
pansion. We take the density profile to decrease as p cx v~ 3 
between v w 0.05 — 0.2 c, as motivated by the typical ve- 
locity v ~ 0.1 c of the dynamically-ejected and wind-driven 
neutron-rich ejecta from NS-NS/NS-BH mergers ( gUJl . We 
assume that p(v) decreases exponentially with v outside this 
range and have verified that our results are not sensitive to 
precisely how we taper the edges of the velocity distribution. 
We find that our results are also similar if we instead assume 
p oc v~ 2 or p cx v~ 4 , implying that the exact density profile 
is not crucial for the overall light curve shape. 

Another input to our calculation is the composition, 
the ionization energies of each element, and the bound- 
bound and bound- free opacities of each element. Unfortu- 
nately, the spectral line information for these very high-Z 
elements is v ery limited. Most of the data available is exper- 
imen tal (e.g. lLawler et aHl2006l 120071 . 120091 ; iBiemont et~ai1 
l2007h . since many body quantum mechanical calculations 
of these elements' spectra represent a formidable task, even 
with modern computing. Much of the experimental work has 
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Figure 4. Bolometric light curve of the radioactively-powered 
transients from NS-NS/NS-BH mergers, calculated assuming a to- 
tal ejecta mass Mtot = 10 — 2 Mq with electron fraction Y e = 0.1 
and mean outflow speed v ~ 0.1 c, and for two values of the 
thcrmalization efficiency ( i|3.2ll , etherm = 1 {solid line) and 
etherm = 0.5 (dotted line). Also shown for comparison (dashed 
line) is a one -zone calculati on b ased on the LP98 model (as im - 
plemented in Kulkarni 2005 and lMetzger. Piro fe Quataertll200gf) 
assuming / = 3 X 10~ 6 (see i|3.2,4H and the same values for Mtot 
and v. 

focused on aiding studies of r-process abundances in ultra 
metal-poor halo stars, which generally make u se of resonant 
absorption lines at optical wavelengths (e.g. ICavrell 1 19961 : 
ISneden et alj|200"3 ). However, the total opacity of most rel- 
evance to merger transients results from densely-packed UV 
lines, for which there is cu rrently insufficient in formation in 
either the Kurucz line list (|Kurucz fc Bel]||l995t ) or the more 
recent experimental studies. Nevertheless, the spectra of at 
least some high-Z r-process elements are almost certainly as 
complex as Fe peak elements, if not more (G. Wahlgren, pri- 
vate communication); this is important because the Fe peak 
elements cause UV "line-blanketing" in normal SN spectra. 
We expect that the same effect is likely to be produced by 
third r-process peak elements since they are largely transi- 
tion metals. 

Given this lack of spectral information, we attempt to 
crudely account for the effects of the unknown r-process 
element lines on the opacity by using the bound-bound lines 
of Fe, but modified to include the correct ionization energies 
of the r-process elements. Specifically, our calculations use 
the ionization energies of Pb as a representative r-process 
element. These uncertainties in the bound-bound transitions 
obviously limit our ability to make detailed spectroscopic 
predictions, but it does allow us to qualitatively address the 
effects of line blanketing on the transients' color evolution. In 
addition, the overall lightcurve shape, the peak luminosity, 
and the characteristic timescale of the event (~ day) are 
robust in spite of these uncertainties. 

4.2 Results 

Figure [4] shows our results for the bolometric light curve for 
a fiducial model with Y e = 0.1, ejecta mass M cj = 10" 2 Mq, 
and outflow speed v = 0.1 c. We show two calculations per- 
formed using different values for the assumed thermalization 
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Figure 5. Top Panel: vh v color light curves from the ethcrm = 1 
calculation in Figure [4] V—, U— , and R— band light curves are 
shown with solid, dotted, and dashed lines, respectively. Bottom 
Panel: Analogous color evolution predicted by the LP98 black- 
body model. 



efficiency, e t horm = 0.5 and ethcrm = 1, which roughly bracket 
our uncertainty in the 7— ray escape fraction ( t|3.2[> . Also 
sho wn for comparis on in Figure [4] is the toy model of LP98 
(cf. lKulkarnJl2005t 323]) , calculated assuming an electron 
scattering opacity and an "f" value = 3 x 10 -6 , as calibrated 
to match the radioactive heating rate in i]3.2.4l Figure [4] 
shows that the light curves predicted using the toy model 
and our more detailed calculation are relatively similar near 
the time of peak emission (f p0 ak ~ 1 day), but their differ- 
ences become more pronounced at earlier and later times. 
The "bumps" in the light curve at t ~ few days in our cal- 
culation are due to recombination of the outer shell electrons 
in our representative high-Z element Pb (and the resulting 
opacity change) as the expanding photosphere cools. 

In the top panel of Figure [S] we decompose the light 
curve into luminosities vL v in several standard optical bands 
(i.e. "colors"). The bottom panel of Figure [5] shows the anal- 
ogous color evolution predicted with the LP98 model, which 
assumes a perfect (single temperature) blackbody spectrum. 
Both calculations predict that the light curve peaks earlier in 
time at shorter wavelengths because the photospheric tem- 
perature decreases with time as the ejecta expands. How- 
ever, the LP98 model predicts an overall vL v peak in the 
UV, while our calculation predicts an earlier peak at longer 
wavelengths (i.e. in V-band) and a clear suppression in UV 
emission at times t J> 1 day. This behavior results from 
strong UV absorption due to dense bound-bound transitions 
( "line blanketing" ) , which produces a much redder spectrum 
than would be predicted by assuming a grey opacity. Indeed, 
rapid reddening following the peak emission epoch is likely 
a defining characteristic of kilonovae. 

We have also explored the sensitivity of our results to 
the mass of the ejecta by performing an otherwise identical 
calculation, but with a lower ejecta mass M ej = 1O" 3 M . 
Our results for the color evolution are shown in Figure [6] 
Although the qualitative features of the light curve evolu- 
tion are similar to the M e j = 10 -2 Mq case, the V band 
light curve peaks somewhat earlier and at a lower lumi- 
nosity, as expected from the analytic scaling relationships 
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Figure 6. Same as the top panel in Figure [5] but calculated for 
M ei = 1O- 3 M . 



ipeak OC M cj 



1/2 



[3]) and Lpcak oc M ( 



1/2 



g]). The 



higher photosphere temperature at the epoch of peak emis- 
sion for lower M c j (r pca k oc M~ 1//8 ; eq. [5]) also results re- 
sults in somewhat bluer peak emission. 



5 DETECTION PROSPECTS 

Because the radioactively-powered emission from NS- 
NS/NS-BH mergers is relatively isotropic, it can in principle 
be detected in at least three independent ways: (1) coinci- 
dent with a source of detected GWs; (2) coincident with 
a short-duration GRB, and (3) via blind transient surveys 
(e.g. PTF and LSST). In this section, we discuss each pos- 
sibility in turn. 



5.1 Gravitational Wave- Triggered Follow-Up 

Advanced LIGO is expected to be sensitive to NS-NS and 
NS-BH merg ers out to a distance ~ 300 and ~ 650 Mpc, 
respectively l|Cutler fc Thorne|[2002]l . For an ejecta mass of 
M c j = 10 _2 Mq we predict a peak V-band luminosity of 
~ 3 x 10 41 ergs s _1 (Figure 3]), corresponding to an absolute 
magnitude My ~ —15. Thus, the entire Advanced LIGO 
volume could be probed by searching down to magnitude 
V ~ 22 — 24. Although this represents a realistic depth for 
a moderately large telescope, the positional uncertainty of 
LIGO/Virgo detections i s expected to ra nge from many ar- 
cminutes to degrees fe.g. lSvlvestrell2003h . As a result, both 
sensitivity and a large field of view (i.e. a high "etendue" ) 
are requirements for any follow-up instrument. Since these 
figures of merit are already optimized for transient survey 
telescopes, projects such as PTF and (eventually) LSST and 
SASIR may be optimal for GW follow-up (in addition to 
their role in blind transient searches; H5.3p . 

Given the short-lived duration ~ 1 day of the ex- 
pected transient signal, rapid GW data reduction (e.g. 
iMarka et ail l2002i ) and dissemination of candidate detec- 
tions to the astronomical community will be crucial for de- 
tection and foll ow-up (as has been adopted by the neutrino 
community; e.g. iKowalski fe Mohrll20Q7l ; lAbbasi et all 2009 ; 



IStamatikos et al.ll2009l ). Indeed, given the unique optical sig- 
nature of NS-NS/NS-BH merger transients 0, moderately 
deep optical/NIR follow-up of even low-significance poten- 
tial GW sources could improve the effective sensitivity of 
Advanced LIGO/Virgo w ith only a relatively mo derate in- 
vestment of resources fsee lKowalski fc Mohrlf2007l for an ex- 
ample applied to high energy neutrino point sources). This 
could prove to be particularly important if the merger rates 
are at the low end of current estimates. Recently, efforts have 
begun to set up a rapid GW data analysis pipeline to pro- 
duce location-probability sky maps within ~ 5 — 10 minutes 
follow ing a GW detection with LIGO/Virgo (|Kanner et al.l 
2008). In fact, a pilot program for the prompt follow-up 
of GW triggers with wide-field optical telescopes such as 
QUEST and TAROT is already underway (P. Shawhan, pri- 
vate communication). 



5.2 Short-Duration GRB Follow-Up 

Given the possible association between short-duration GRBs 
and NS-NS/NS-BH mergers, another method to search for 
kilonovae is with deep optical/IR observ ations on ~ 1 
day ti mescale following the burst. Recently, iKocevski et al.l 
(2009) presented upper limits on the presence of a putative 
LP98 kilonova using follow-up observat ions of GRB 070724A 
and GRB 050509b l|Hiorth et al.ll2005l ). Our results in Sf4J2] 
show that the LP98 model (properly calibrated) does a rea- 
sonably good job of reproducing the qualitative features of 
the optical/NIR light curves around the time of peak emis- 
sion. Thus, by assuming v w 0.1 c and / w 3 x 10" 6 f i)3.2.4|l 
we conclude from their Figures 8 and 9 that M e j must be 
<; 0.1M & and £ 10~ 3 M Q for GRBs 070724A and 050509b, 
respectively. Although the former (070724A) is not partic- 
ularly constraining on merger models because such a large 
ejecta mass is unlikely, the latter non-detection (050509b) 
suggests that the disk that formed in this event was rather 
small (< 10~ 2 M Q ; see the discussion in A low disk 

mass is, however, still consistent with a merger interpreta- 
tion for this event because the isotropic lumin osity of the 
GRB was quite low (£ 7l i so ~ 2.4 x 10 48 ergs; iKann et all 
2008): even ignoring geometric beaming corrections, only 
an accreted mass ~ 10 _5 Mq is necessary to produce a rel- 
ativistic jet with energy i? 7 ,i so assuming that the efficiency 
for conve rting rest-mass e nergy to 7— ray power is ~ 10 per 
cent (e.g. lMcKinneyll2005l ) . 

iBerger et al" l (|2009l ) present additional early opti- 
cal/NIR follow-up observations of GRB 070724a, including 
the discovery of transient emission peaking t ~ 3 hours fol- 
lowing the burst. Due to the transient's very red spectrum 
(which is highly unusual for a standard GRB afterglow) 
they suggest a possible NS merger transient interpretation. 
They conclude, however, that this possibility is unlikely: the 
brightness and rapid rise time of the transient require val- 



ues of / ~ 5 x 10~ 3 and M ei ~ 10 Mq which, within the 
standard LP98 model, predicts a photospheric temperature 
peaked in the UV (eq. [3]) and thus inconsistent with the 
transient's red colors. Our calculations in ^4]suggest that UV 
line blanketing could produce a redder spectrum (see Fig. [5}. 
However, the value of / ~ 5 x 10 -3 they nominally require 
is three orders of magnitude higher than we predict from 
radioactive decay ( i]3.2.4[l . It thus appears more likely that 
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the emission detected by Berger et al. is afterglow-related 
(with the observed reddening perhaps due to dust). 

Perhaps the most promising candidate kilonova detec- 
tion to date was following GRB 080503, which showed an 
unusual rise in its optical afterglow light curve at t ~ 1 day 
befor e rapidly fading over the next several days |Perlev et al.l 
2009). Although limited color information was obtained near 
the emission peak, the observed light curve evolution is 
largely consistent with that expected from a kilonova for 
an assumed ejecta velocity v ~ 0.1 c and mass M e j ~ few 
x10~ 2 Mq. However, although the event was well-localized 
on the sky, no obvious host galaxy was detected coincident 
with the burst, despite the fact that a relatively low red- 
shift z ~ 0.1 would be required to fit the observed peak 
brightness (if the event was indeed a kilonova). In princi- 
ple, NS-NS/NS-BH binaries can receive natal "kicks" from 
their supernovae which may eject them into intergalactic 
space, thereby making it difficult to identify their original 
host galaxy. This possibility is consistent with the very low 
density of the circumburst medium inferred for 080503. Per- 
haps more problematic for the kilonova interpretation, how- 
ever, is the X-ray detection by Chandra coincident with the 
optical rise, which suggests that the optical emission at ~ 1 
day may simply be due to an (albeit unusual) non-thermal 
afterglow. The fact that the X-ray luminosity exceeds the op- 
tical luminosity appears especially difficult to explain if both 
are related to radioactively powered quasi-thermal emission. 

Although the possible presence of kilonovae following 
short GRBs is not well-constrained at present, this situ- 
ation could in principle improve with additional sensitive, 
early-time optical/NIR observations of short GRBs. Indeed, 
Swift should remain operational through the next decade 
and presently detects short GRBs at a rate of about ~ 10 
per year. Unfortunately, however, this approach may en- 
counter fundamental obstacles due to the non-thermal af- 
terglow emission which also generally accompanies GRBs. 
In most accretion-powered GRB models, the lu minosity of 
the burst increases with the accretio n rate (e.g. McKinnevl 
120051 ; [ Zalamea fe Beloborodovl [2009) and, hence, with the 
disk mass. Since the quantity of neutron-rich ejecta may be 
a relatively constant fraction of the disk mass (|Metzger et al.l 
2009a; H2.ip . the luminosity of the kilonova (L pea k oc M^/ 2 ; 
eq. [4[) may positively correlate with the luminosity of the 
GRB(j Since the afterglow luminosities of short GRBs ap- 
pears to scale with the prompt GRB fluence (as in long- 
duration GRBs; [Nvsewander et all 2009), the afterglow may 
generically swamp any putative kilonova emission. To date, 
this appears to be true even in cases when the circumburst 
densi ty appears to be ve ry low and the afterglow is relatively 
dim l|Perlev et al.ll2009l) . 



5.3 Blind Optical Transient Surveys 

In this section we assess the prospects for detecting kilono- 
vae from NS-NS /NS-BH mergers with present and upcoming 
optical transient surveys. The virtue of such a search strat- 



3 An exception may occur in the case of BH-NS mergers, where 
in some cases large amounts of material ca n be ejected rel ative to 
the mass of the accretion disk that forms (Rosswog 2005) 



egy is that it does not rely on a GW or high-energy EM 
trigger. 

Based on observed binary NS systems, iKalogera et al.l 

( 2004) find that the NS-NS merger rate in the Milky Way is 
between 1.7 x 10~ 5 and 2.9 x 10" 4 yr" 1 at 95% confidence 
Population synthesis estimates (e.g.. lBelczvnski et al.1120061 ) 
are consistent with this range but with larger uncertain- 
ties. Since there are no known BH-NS binaries , the B H-NS 
merger rate is even less certain. [Bet he fc Brownl l|l998T ) argue 
that BH-NS mergers coul d be substantially m ore common 
than NS-NS mergers, with I Bet he et all (|2007l ) estimating a 
rate ~ 10 4 Gpc _3 yr _1 , corresponding to ~ 10 -3 yr -1 in the 
Milky Way. An interesting limit can be placed on the total 
amount of neutron-rich ejecta from NS-NS/NS-BH merge r 
from Galactic chemical evolution fe.g. lMetzger et aTll2009al ). 
Accounting for the total observed abundances of elements 
with A > 100 in our Galaxy, for example, requires an average 
production rate ~ 10~ 6 M Q yr" 1 (e.g. lQianll2000l ). Assum- 
ing that a merger ejects M e j ~ 10 _2 (10 _3 )Mq on average, 
the Galactic merger rate cannot exceed ~ 10~ 4 (10~ 3 ) yr -1 
in order to avoid over-producing these rare neutron-rich iso- 
topes (see further discussion in H5.3.1l below). 

Assuming that the NS-NS rate is p r oportional to the 
blue stellar luminosity (|Phinnevl Il99ll ; iKopparapu et al.l 
2008), a Galactic rate of .Rns-ns = 10 _4 7?_4 yr -1 corre- 
sponds to a volumetric rate of 10 _6 i?_4 Mpc -3 yr -1 . For 
v — 0.1 c and M e j ~ 10 _2 M e j,_2MQ our calculations predict 
an optical transient with a peak V-band luminosity vL v ~ 
3 x 10 41 Af 1 / 2 2 ergs s _1 (Fig.[S]and[S|. For a limiting magni- 
tude Mv = 25(24) [21], this corresponds to a maximum de- 
tection (luminosity) distance of D L = 1070(680) [170]MV1 2 
Mpc and a co-moving volume V = 2.8(0.9) [2 x 10~ 2 ]M 3 / 4 _ 2 
Gpc 3 . The Palomar Transie nt Factory (PTF) 5-day ca- 
dence survey (|Law et al.ll2009l ), which surveys an active area 
~ 2700 deg 2 to a limiting AB magnitude of 21, should there- 
fore detect ~ 1.4J?_ 4 M 3 / 4 . 2 yr" 1 . Thus, if NS-NS mergers 
occur at the upper end of present rate estimates (~ 10~ 4 
yr -1 ) and M e j ~ 1O _2 M0 is indeed representative, current 
surveys such as PTF should "blindly" detect ~ 1 merger per 
year. We emphasize, however, that the total predicted rate 
of events and their luminosity function is sensitive to the 
distribution of ejecta masses M e j, which in principle could 
range from ~ — 0.1 Mq given present uncertainties ( i)2,ip . 

Prospects for detection are much better with the Large 
Synoptic Survey Telescope (LSST), which will image the 
entire sky down to a limiting magnitude ~ 24.5 every 3-4 
nights and should detect NS-NS merger events at a rate 
~ 2 x 10 3 7?_ 4 M 3 / 4 2 yr -1 . Note that LSST is expected 
to come on-line in 2015, roughly coincident with Advanced 
LIGO/ Virgo, and together they have the potential to com- 
pletely revolutionize our understanding of compact object 
mergers. 

Other thermal transients are predicted to occur in Na- 
ture on ~ day timescales, which could be confused with 
NS-NS/NS-BH mergers. Examples include ".la" SNe due to 
unst able thermonuclear He flashes from white dwarf bina- 
ries (|Bildsten et al.ll2007i ; IPoznanski et al.ll2010l ) and Nickel- 
rich outflows from th e accretion-induced collapse of WDs 
l|Metzger et al.ll2009bl ; Darbha et al. 2010). Such events may 
originate from a similar stellar population to NS-NS/NS-BH 
mergers. However, one "smoking gun" feature of kilonovae 
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Figure 7. Rates of detected kilonovae from NS-NS/NS-BH merg- 
ers with present and upcoming surveys as a function of the merger 
rate iV me rge (bottom axis) or the average ejecta mass M e j (top 
axis), calculated under the assumption that NS-NS/NS-BH merg- 
ers are the primary source of r-process elements in the Galaxy. 
Also shown are the NS-NS merger rate estimates (95% confidence 
interval) from Kalogera et al. (2004). 



from NS-NS/NS-BH mergers is the presence of optical ab- 
sorption lines due to heavy neu tron-rich elements ( as in re- 
process enriched halo stars; e.g. ISneden et al ] l2003h . which 
may not be present in white dwarf systems (although some re- 
process nuclei may be produced in neutr ino-heated winds in 
the case of AIC; e.g. lDessart et alj|2006h . Thus, NS-NS/NS- 
BH mergers may be distinguishable from other transients 
with rapid, deep spectroscopic observations. We plan to ex- 
plore more detailed calculations of these r-process spectral 
features in future work. 

In reality, only limited information will initially be avail- 
able to transient searches (e.g. photometric colors, at best). 
Thus, the reddening of kilonovae at times > t pea k (Fig.[SJ), if 
indeed robust, may be crucial for identifying these events. In 
fact, of the variable s ources from t he Slo an Digital Sky Sur- 
vey characterized by ISesar et al.l (|2007l ). only a small frac- 
tion are as red as we predict kilonovae from NS-NS/NS-BH 
mergers to be at peak light (U— V > 2 following p eak bright- 
ness: compare our results in Fig. [5] with Fig. 4 of lSesar et al.l 
2007). This suggests that a promising search strategy for 
detecting kilonovae is to trigger on anomalously red events 
of duration ~ hour— day for more detailed follow-up obser- 
vations. 



5.3.1 Implications for the Origin of r -Process Elements 

The astrophysical origin of the r-process elements re- 
mains one of the great mysteries in nuclear astro- 
physics (see lQian fc Wasserburell2007l for a recent review), 
with the t wo chief candidat es being core-collapse super- 
novae (e.g. iMever et al"1ll992T ) and NS-NS/NS-BH mergers 
(|Freiburghaus et al. 1 1 19991 ; see Fig.©. Since the luminosity 
of merger transients is directly related to their nucleosyn- 
thetic yield (eq. [3]), this implies that the detection of, or 
constraints on the rate of, kilonovae from NS mergers di- 
rectly probes the origin of r-process elements. 

As a concrete example, if one assumes that NS mergers 



are the dominant source of r-process elements in our Galaxy, 
then the mean mass ejected per event M e j and merger rate 
iVmerge are related by iV mcrg e = 10~ 4 yr- 1 (M cj /lO _2 M )- 1 , 
where we have assumed a Galactic r-process production 
rate of 1O" 6 M yr" 1 (e.g. Qian 2000). Since the kilonova 

1/2 

luminosity L pca k oc M e - , the rate of detected transients 
oc L 3 l 2 N merge oc M c j 1,/4 oc iVmorg C . Figure [7] shows the ex- 
pected detection rates versus Af me rge (or, equivalently, M e j) 
with present and upcoming transient surveys if one assumes 
that NS-NS/NS-BH mergers are the dominant r-process 
source. Note that within the current uncertainties in Af me rge, 
current transient surveys should detect a few events per year 
if mergers are indeed the dominant source of the r-process, 
independent of M e j. In reality, the detection efficiency for 
low Mej may be somewhat lower than this simple estimate 
due to the shorter transient duration £ pea k oc MV 2 (eq. [3])- 



6 CONCLUSIONS 



In the ir seminal paper on nucleosynthesis, iBurbidge et al.l 
l| 19571 . B 2 FH) proposed th at SNe are powered by the ra- 
dioactive decay of 254 Cf (cf. IBurbidge et al.lll956t ). Although 
we now appreciate that most supernovae are powered by 
56 Ni and 56 Co, the B 2 FH picture of "r-process-powered" 
SNe still holds relevance for the neutron-rich ejecta from 
NS-NS/NS-BH mergers. 

In this paper we have presented the first calculations 
of the radioactively-powered transients from NS-NS/NS- 
BH mergers that self-consistently determine the radioactive 
heating using a nuclear reaction network and which accu- 
rately model the light curve and color evolution with a ra- 
diative transfer calculation. Our main conclusions are sum- 
marized as follows: 

• The radioactive heating E on timescales t pca k ~ 
hours— days results [in approximately equal parts] from the 
fission and ft— decays of heavy nuclei, which are produced 
by the r-process at much earlier times (t <, 1 second; see 

Fig- HI). Our results for E at t ~ t pca k are relatively insensi- 
tive to the precise electron fraction and early-time expansion 
of the ejecta (e.g. whether it is dynamically-ejected or wind- 
driven), and to details of uncertain nuclear physics such as 
the theoretical nuclear mass model. 

• The net heating rate decreases approximately as a 
power law Q oc t~ a with a ~ 1.1 — 1.4 for t ~ hours— days, 
similar to the assumption Q oc t" 1 in the LP98 model. The 
total heating rate is ~ 3 x 10 10 ergs s _1 g _1 at t ~ 1 day. 
By calibrating the LP98 model using our results, we find an 
effective "f" parameter ~ 3 x 10 -6 which is generally much 
lower than previously assumed. 

• fi— decay electrons ( H3.2.1[l and fission daughter nu- 
clei ( i]3.2,2[) both thermalize with the plasma on timescales 
~ £ pC ak, while only a portion of the 7— rays likely thermal- 
ize. We estimate that the net thermalization efficiency is 

etherm -0.5-1 (MM ■ 



For an ejecta mass M e - 



10~ 2 A/ Q to 10 _d M, 



predict a transient that peaks on a timescale 
a bolometric and V-band luminosity - 1 " 1 ~ 
= -16) and ~ 3 x 10 41 ergs s _1 (M v 
(Figs. H and [g. 



we 
1 day at 
10 42 ergs s" 1 (M B 
-15), respectively 
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• We argue that the transition metal r-process elements 
are likely to have UV absorption due to line blanketing (like 
Fe peak nuclei). As a result, we predict that NS merger 
transients will be relatively red (and redden in time; see 
Figs.[5]and[6]), a prediction not captured by assuming single- 
temperature blackbody emission. More detailed models of 
the color evolution of NS merger transients will require a 
better understanding of the UV and IR spectral lines of sec- 
ond and third r-process peak elements ( H4.ip . The presence 
of absorption lines due to heavy r-process elements is one 
"smoking gun" prediction of NS merger transients. 

• Because NS merger transients are isotropic, they can in 
principle be detected in three independent ways: in coinci- 
dence with a detected GW source; following a short-duration 
GRB; and with blind optical/NIR transient surveys. 

• Given the low luminosities and rapid evolution of kilo- 
novae from NS mergers, their detection will require close 
collaboration between the GW and astronomical commu- 
nities. Given the unique observational signature of kilono- 
vae, the real-time follow-up of GW detections with sensitive, 
wide-field telescopes could improve the effective sensitivity 
of LIGO/Virgo. 

• For an average ejecta mass M c j m 10~ 2 Mq, current 
surveys such as PTF should "blindly" detect ~ 1 NS merger 
transient per year if the merger rates lies at the high end of 
present estimates; LSST should detect ~ one thousand per 
year under the same assumptions. 

• Since the luminosity and detection rate of NS merger 
transients is closely related to the yield of heavy neutron- 
rich elements, current transient surveys are directly probing 
the unknown astrophysical origin of the r-process (Fig. 0). 
Holding the total r-process injection rate in the Milky Way 
fixed at 10 _6 Mq yr _1 implies a detection rate ~ few yr _1 
and ~ 10 3 yr _1 for PTF and LSST, respectively, indepen- 
dent of the average ejecta mass. 
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